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Abstract 

Jahn- Teller (JT) effect both in the absence and presence of the strong 
Coulomb correlation is studied theoretically focusing on the reduction AK 
of the kinetic energy gain which is directly related to the spin wave stiff- 
ness. Without the Coulomb interaction, the perturbative analysis gives 
AK I (g 2 1 MQ, 2 ) = 0.05 — 0.13 depending on the electron number [g: electron- 
phonon(el-ph) coupling constant, M: mass of the oxygen atom, fi: frequency 
of the phonon]. Although there occurs many channels of the JT el-ph inter- 
action in the multiband system, the final results of AK roughly scales with 
the density of states at the Fermi energy. In the limit of strong electron cor- 
relation, the magnitude of the orbital polarization saturate and the relevant 
degrees of freedom are the direction (phase) of it. An effective action is de- 
rived for the phase variable including the effect of the JT interaction. In this 
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limit, JT interaction is enhanced compared with the non-interacting case, and 

AK is given by the lattice relaxation energy El for the localized electrons, 

although the electrons remains itinerant. Discussion on experiments are given 

based on these theoretical results. 
71.27.+a, 75.30.-m, 75.30.Et 
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I. INTRODUCTION 



It has been recognized that in the strongly correlated electronic systems both the 
electron-phonon (el-ph) and electron-electron (el-el) interaction are enhanced and play im- 
portant roles. In manganites, the colossal magneto-resistance has been discussed 
from both points of view. In this system the ferromagnetism is basically explained by the 
double exchange model. However the orbital degeneracy of the e g -states is considered to 
be very important, and the orbital ordering or disordering is the crucial issue for the under- 
standing of the CMR effect. Therefore there are two viewpoints on this problem, namely 
the orbital degrees of freedom is governed by (i) the Jahn- Teller (JT) electron-phonon cou- 
pling or (ii) the electron correlation. In the former case, the change of the bandwidth due to 
the crossover from small to large JT-polaron is the key mechanism of the CMR effect, and 
the JT effect is assumed to be negligible in the ferromagnetic metallic state.@~il This pic- 
ture appear contradicting with the orbital orderings, which require rather strong coupling, 
surrounding the ferromagnetic region in the phase diagram.01lhH However it might be the 
case that the metallic screening weakens the el-ph interaction and/or the el-el interaction, 
and only in the ferromagnetic metallic state both of them could be neglected although the 
Hund's coupling is strong enough to polarize the spins perfectly. On the other hand, it has 
been recognized by several authorsiHilIzlIl that the Coulomb interaction play the important 
role in the physics of the orbital degrees of freedom, and the JT-interaction is the secondary 
effect. 

In this paper we revisited this issue by considering both the el-ph and el-el interaction, 
because both of them are considered to be relevant. Then the interplay between these 
two interactions is the key issue. As shown below, the metallic screening of the JT el-ph 
interaction does not occur in contrast to the coupling between the breathing mode and the 
charge fluctuation. Coulomb interaction and JT effect collaborate with each other, and 
it is concluded that JT effect is enhanced by the el-el interaction by comparing the two 
limits of zero and strong electron correlations. Similar idea has been proposed by one of 
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the authorsa in the context of the large-c? approximationEJt2rEJ, where the single site model 
embedded in the dymanical environment is considered!!. To our knowledge the effect of the 




JT phonon scattering with d = 3 has not yet been exploited so much, and we found that 
the el-ph coupling which modulates the transfer integral becomes relevant in both limits 
of the non-interacting and the strong correlation limit. Another important issue is how 
the el-ph intereaction manifests itself in the spin stiffness observed in neutron scattering 



kinetic energy.f-3 The polaronic effect on the kinetic energy is therefore expected to appear 
as a reduction of the stiffness.0'il. The apparent absence of this reduction has led to the 
conclusion that el-ph interaction is irrelevant in the ferromagnetic metallic state. However 
we found that JT el-ph interaction is there even in the ferromagnetic state although it might 
be hidden in the inaccuracy of the theoretical estimation of the bare spin stiffness compared 
with the experimental ones. We expect near 3% of the reduction of the kinetic energy and 
spin stiffness due to the JT interaction as a lower bound estimated in the non-interacting 
limit. 

We first consider the non-interacting electrons with orbital degeneracy with JT el-ph 
coupling.0 The single-band model where the charge density is coupled to the phonon can 
be treated by means of the canonical transformation, leading to the reduction of the kinetic 
energy via the Debye- Waller factor.@ This argument is applicable to the interaction with 
the breathing mode. However this method can not be generalized to the JT el-ph interaction 
with multiband electronic structure because it includes off-diagonal components with respect 
to orbital indices. Due to this difficulty, it is hard to apply the same argument as the 
breathing mode case to clarify whether the JT polaron also leads to the reduction of the 
kinetic energy or not. Therefore we employ the perturbative analysis on the JT el-ph 
coupling to estimate the reduction of the kinetic energy and spin stiffness. Although many 
channels contribute, each of which can be even negative, the resultant reduction in the 
kinetic energy gain is roughly proportional to the density of states at the Fermi energy, and 
has a peak at around n = 0.6 and n = 1.4, where n denotes the filling of the e g band. 




is roughly proportional to the 
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Next we consider the strong correlation limit by employing the effective Lagrangian which 
is derived as a projection onto the polarized orbital state.Hl In this limit, the magnitude of 
the orbital polarization has been saturated, and the only degrees of freedom coupled to the 
JT phonon is therefore the direction of the polarization (which corresponds to the shape of 
the orbital ordering). The direction is determined by the double exchange interaction, which 
is of the order of the transfer integral, to- The quantum fluctuation of the direction (physi- 
cally the fluctuation of the orbital shape) is couples to the lattice deformation via the el-ph 
interaction, leading to the lattice relaxation energy El = g 2 /MQ 2 (g, M and Q denote the 
JT el-ph coupling constant, the atomic mass, and the phonon frequency, respectively). The 
characteristic frequency u n of the orbital fluctuation is of the order of the transfer integral 
to, and is larger than the phonon frequency Q. Namely the orbital deformation can follow up 
the phonon, leading to the kinetic energy correction, AK ~ El- The phonon frequency 
dependence of the kinetic energy correction in this case, [AKu-mo /-Ex] ~ O (1), is differ- 
ent from that of non-interacting electrons case, [AK U=0 (fl) / E L ] ~ (fi/to) / (1 + Q/t ) 2 . 
Thus the strong correlation enhances the JT effect, in sharp contrast to the case of breath- 
ing mode where the Coulomb interaction reduces the el-ph interaction. This is understood 
rather easily. In the case of the breathing mode, the Coulomb interaction suppresses the 
charge fluctuation while the breathing mode induces it, i.e., these two interactions compete 
with each other, and the former suppresses the latter. Furthermore the metallic screening 
effect also suppresses the charge fluctuation (This situation has been discussed in the context 
of the vertex correction of the el-ph interaction in the physics of high-T c cuprate M). On the 
other hand, in the JT mode case, where el-el and el-ph interactions collaborate to induce 
the orbital pseudospin moment, the former enhances the latter and vice versa. Therefore 
there is no reason to expect the weakening of the JT el-ph interaction with the doping when 
the strong el-el interaction keeps the orbital pseudo-spin moment to saturate even in the 
ferromagneitc metallic state. 

The plan of this paper follows. The perturbative analysis of the JT el-ph interaction 
for the noninteracting electrons with orbital degeneracy is given in Section II. The strong 



correlation limit is studied in Sec. Ill, and discussion and conclusions are presented in Sec. 
IV. 



II. NON-INTERACTING LIMIT 

The JT interaction is given as 

Hjt = 9^2 [(d-ladja — dj b dj b ^ ■ Q u j + (d^ ja djb + d^ b dj^) ■ Q v j , (1) 

j 

with a coupling constant g. The spinless operator for the half-metallic ferromagnetic phase, 
<ij 7 , creates a spin polarized e g electron with orbital 7 [= a(d x 2_ y 2), b(d 3z 2_ r 2)] at site j. 
Q u and Q v denote the normal coordinates of the displacement of the oxygen ions A a (a = 
x,y,z): Q u = (2A Z — A x — A y ) /v6, Q v = (A x — A y )/^/2. Let us consider the kinetic 
energy correction due to the JT phonon scattering with a two-band model, 



i(5,77' 
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1 - - Mil 
2M Pj ' Pj + ~lT Ql ' Qj 



+ H JT . (2) 



are realistic anisotropic hopping intensities given in ref.0. Qj is defined as 
(Quj-, QvjY- M, f2, and Pj denote the atomic mass, the elastic constant, and the canonical 
momentum of Qj, respectively. For a simplified model without orbital indices, the canonical 
transformation is a standard method to deal with the electron-phonon interaction, leading 
to the kinetic energy reduction by the Debye- Waller factor.^! The canonical transformation 
for the JT phonon has off-diagonal elements with respect to orbital indices, due to which 
the application of this method becomes complicated. We therefore employ the perturbative 
calculation of the kinetic energy K up to the second order with respect to the coupling 
constant g, as 

K = K + AK , AK ~ O (g 2 ) . (3) 
AK is expressed by diagrams of the self-energy shown in Fig. [I]. 
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FIG. 1. Diagrams of the self-energy due to the JT interaction. Correponding to the doubly 
degenerate orbitals, there are two kinds of vertices and the diagram is composed with four-channeled 
propagator. 



There are two vertices corresponding to Q u - and (^-scatterings. The propagator takes 2x2- 
matrix form with respect to the orbital indices. Annihilation and creation operators of the 
JT phonons are introduced as 



Qu, v (t) = J ^ ■ (a u , v ■ e~ int + al v ■ e**) . 



. 2Mtt 

With the momentum representation of the operators c = (a, d), 

1 



(4) 



Cj (r) 



iuji ->• z , 



(5) 



with Matsubara frequency iu>i, the propagators of electrons and phonons are given as 



9? (z) = -T (d kl (z) d[ Y (z)) o itj z = [{z + n) 5 W - e? 
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T (a qMv) (z) a\ Mv) (z))+T (al qMv) (z) a_ qMv) (z)) 
1 ' 1 =D q (z) , 



z + f2 z — Q 



(7) 



respectively. The u- and -u-modes have the same mass and frequency because they belong 
to the same irreducible representation, and then the same phonon propagator. Here we 



neglected the inter-cluster coupling of the vibration and hence the wavelength dependence 
of the phonon frequency. Coefficients A±, k in Eq. (|j) are defined as 
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(8) 



with dispersion relations of the hybridized bands given as 
H (±) _ 1 

—ft — n 



(9) 



The orbital index 7 is used as a = 6 and b = a. ef 1 (= eg), e| 6 e|) , and eg 6 are the cosine 
dispersions with overlap integrals between the orbitals \a) = \x 2 — y 2 ) and \b) = \3z 2 — r 2 ). 
With the propagators, the kinetic energy correction AK is given as 
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^ K = o E E E 4 7 <?f 1 • £? 72 • 9j f • [ujo , 

i ft,77' 7172 



(10) 
(11) 



where the contour c surrounds the poles of the fermi distribution function / (z). In Eq. (11), 
g^t (z) and g^Jg (z) correspond to the scattering by u- and vert ex, respectively. These 
contributions are represented by diagrams shown in Fig. ^. 





FIG. 2. Corresponding diagrams of the terms in Eq. (11). 



Introducing a notation, 



<*) = e s r : <*■*> • s? 7i i**) , (12) 



Eq. fllOD is expanded as, 
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+1 V [e 66 • / (66) 

^ Z^[ E k ' 1 (ab;aa) £ fc ' J (a6;6fe) £ fc ' ^(ab;afe) 1 {aa;bb) )\ ' \ LO ) 
k 

Three terms correspond to the contribution from ££ a , E| 6 (= S^ a ), and E^ 6 , respectively. 
Note that the upper (lower) suffix 7«7^ of -f^y.^y ^ means that the correpsonding contribu- 
tion comes from a diagram composed of a propagator g^" 7 " (g 7 ™ 7 ™) for the state \k — q) (\k)) 
[When one finds (ab) in the upper (lower) suffix, that contribution contains the hybridiza- 
tion during the propagation with the wave vector \k — q) (\k))\. In Appendix A are given 
concrete forms of // 7171 ^ , v 

Numerical results are shown in Fig. ^| with parameters to = 0.72 evil and Q = 0.05 
eV.El The kinetic reduction AK is calculated as a function of the filling n. 
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FIG. 3. Kinetic energy reduction as a function of the filling n for non-interacting electrons, (a) 
Total result and partial contributions from u- and f-vertices in a realistic two-band model with 
e g anisotropy. (b) Result for one and three dimensional single band cases with breathing type 
phonons. (c) Density of states as a function of the fermi energy for given n (lower n-axis) in the 
system used for (a). 



Panel (a) shows the total result [Eq. (0)] and partial contributions due to the u- and v- 
vertices. The particle-hole symmetry with respect to the axis n = 1.0 is seen. The positive 
definite result is obtained for the whole range of n. In order to understand the origin of the 
n-dependence in the plot (a), we also calculate the simpler case with breathing type phonons 
and single band electrons [Eq. ( |A8|) 1 in one and three dimensions, as shown in the panel 
(b). In this case AK/El roughly scales to the density of states at the fermi level N(ef) 
for given n (see Appendix B). According with this expectation, the result has the minimum 
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(maximum) at n = 0.5 for one (three) dimensional case. The positive definite result in the 
plot (b) is consistent with the consequence from the canonical transformation method.ll 
Though the expression for the realistic doubly-degenerate e g case [Eq. flTJ|) and the panel 
(a)] is much more complicated, the result also seems to scale to the density of states. For the 
comparison, the plot of the density of states in this case is shown in the panel (c). Values of 
n giving the peak and the dip in panel (a) and (c) actually coincide each other. As discussed 
in Appendix B, k- (wave vector) points near the Fermi level contribute dominantly to AK. 
Unless the n-dependence of each contributing value is so sensitive, AK simply scales to the 
population of the contributing fc-points and hence N(ef (n)) for given n. This gives a rough 
explanation for the correlation between AK (n) and N(ep)- More intuitively, AK scales 
to the population of electrons around the fermi surface [oc N (ep)] which is subject to the 
phonon scattering. 

In order to see how each scattering process contributes, we re-divide AK into several 
contributions as 

— y 2-^1 k ' 1 {aa ] aa) + £ k ' 1 (bb;bb)\ + 9 ^ I k ' (a b ;ab) + £ k ' 1 (ab;ab) 
k k 
, \ c ah ( j{aa) j{aa) \1 ^-v \ aa j{ab) , 66 Aab) 

2^ l b k ' V J (a6;aa) ~T J (a6;66) )\ ^ [ b k ' 1 (ab;aa) ~T b k ' 1 (ab;bb) 
k k 

+E[*r('SV / £W] > ( 14 ) 

k 

and plotted each contribution separately in Fig. B. 
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FIG. 4. Filling (n-) dependence of the each contribution defined in Eq. (|j~4|). 



The first, second and third terms are coming from E 7 / 7 , whereas the fourth and fifth from 
Ej 7 . The latter contribution is small compared with the former. The first term, which 
is diagonal with respect to the suffices of /, just corresponds to the superposition of two 
breathing type diagrams with u- and v-vertices, respectively. (Remember the suffix rule of 



the JT interaction, which are sorted further into two classes. One class (the second and 
fourth terms ) corresponds to the twice inversion of the orbital state (correspondingly the 
off-diagonal orbital suffix ab appears twice) to come back to the original orbital state (like 
a — > b — > a). Consequently this class picks up the diagonal dispersion ej? as its weight. The 
other class (the third and the fifth terms) with odd number the orbital index ah thus picks 
up the off-diagonal weight (like a — > a — > b). Because e^ b roughly corresponds to the 
energy scale of the stabilization due to the band hybridization, it gives basically the negative 
contribution (stabilization of the energy) as shown by the behavior of the third term in Fig. 

This stabilization around n — 1 is mainly attributed to the dip with negative values of 
the "3rd-term" seen in Fig. f| (a). This behavior is understood as follows. The third term 
oc e% b reflects the stabilization of the lower band due to the repulsion with the upper band 
(its magnitude is tf b i+s ). Such a stabilization is most remarkable at the region where the 
hybridizing two bands cross with each other. For the half-filled case, n = 1, the fermi level is 
located at the middle of the bandwidth, where the band-crossing occurs, leading to the most 
effective stabilization. That is why the stabilization of the third term is most remarkable 
around n — 1. 




mentioned before). The other terms arise due to the multiband structure and 



III. STRONG CORRELATION LIMIT 



Strong on-site repulsions in e g orbitals can be written asEH 

Hoxi.— site P ^ ^ Tj • Tj , Tj — ^ ^ dj^c^idj^' , 




(15) 



j 
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with spinless operators. Tj is the isospin operator representing the orbital degrees of freedom 
with 2x2 Pauli matrices cx 7 y = (cr^y , cr^y , cr*y J . (3 is a parameter of the electron-electron 
interaction of the order of the Hubbard repulsive U. This interaction induces the finite 
orbital polarization, which can be represented by the Stratonovich-Hubbard field (orbital 
fluctuation field) (p T as0 



H e l = E d h & - I*) d n + E tlj'+sdl^+S^' + E 



ii5,77' 



1 s • .\/o 2 , 



2M 



-<5j ■ + {30t 



(16) 



It is seem that the orbital fluctuation field tp? as well as the JT phonon Q is coupled to 
the isospin T in the form of linear combination 2f3ip = 2(3cpx ~ gQ$- After integrating out 
phonon coordinates, the effective action in terms of the field ip is obtained as 



S, 



eff 



E 4j'7 (^t- — A 4 ) 4/7 + X! tfJ+sd^di+Sa' 

hi i<5,77' 

2/3 ^ f iin • ^> • 



2(3M (qj + Q 2 
2/3M(u; 2 + ft 2 )+2 



2 r ^ 



(17) 



where cj n = 2im/f3 is the Matsubara frequency for the bosons. The phonon dynamics 
induces the retardation effect for the field if), which is represented by the a; n -dependence of 
the second term in the above equation. Now let us assume that the electron correlation is 
much larger than the JT coupling, namely (3 3> El- It is noted here that we do not assume 
El *C ^cb namely the weak coupling limit. Then we can expand in the JT coupling g in 
eq.(17) as 



S, 



eff 



dr 
+/?E 



E 4?7 ~ ^) dj-y + E tJJ+sdiydi+srf 

3,7 j<5,77' 
2 
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2/3M(w 2 +ft 2 )_ 



•^,n'0i,n-2/3E T i,n-'0. 



(18) 



Because we are now interested in the strong correlation limit, (3 ^> t Q , the magnitude of the 
orbital polarization is fully developed. This corresponds to the fixed l^l = Pt = 1/2, and 
we consider its direction only within the xz-p\ane because the JT coupling prefers the real 
orbital states^. Then ijj is parametrized as 
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4j = V • * (sin 0-,-, 0, cos , (19) 

with the phase angle 9j being the only relevant degrees of freedom. Correspondingly, the 
isospin is forced to be parallel to ip, and hence the Grassman variables c?] T ,c?j 7 are replaced 
by 

4 7 = cos(%/2), /t sin(%/2)] , (20) 

with the spin/orbital-less fermion variable p, f. Putting this expression into Eq. fll8f) , the 
kinetic energy term can be written as 

^2t iti+s (di, 9 i+s ) ■ fifi+s , (21) 

i,S 

with the ^-dependent transfer integral, 

ti,i +5 {9 h 9 i+5 ) = t\l +6 cos^/2) coa{6 i+s /2) + tf i+5 sin^/2) sm{9 i+s /2) 

+ t\ 2 t+5 cos(^/2) sin(^ ;+5 /2) + tf t+5 sin^/2) cos(9 i+5 /2) . (22) 

This gives the coupling of ^-field to the fermion. On the other hand, the dynamics of the 
^-field is generated through this coupling by integrating over the fermions f \ f. 

S° cS = n(<f, u n ) ■ 9(q, iu n )9(-q, -u n ) , (23) 

where 9 is measured from the mean field value, and H(q, u n ) is the orbital correlation function 
of the fermions. Although the quantitative results depend on the details of the model and 
the orbital ordering to start with, the orbital fluctuation has a gap of the order of to and there 
occurs no infrared divergence. Therefore the characteristic frequency uj n and wavevector q 
for H(q,u n ) are ~ to, and 7r/a (a: lattice constant), respectively. The Hamiltonian consists 
only of the kinetic energy besides the JT coupling term, and its expectation value is given 
by 

<*•> =W = 0,«w = 0) (/,*/«)„ M + fckk>n(,» . (24) 

i,S " q,u 
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Now let us analyse the correction term due to the JT interaction 



,2 



jr 

jO 2 



ss °s = ~J2 2M (uj 2 + tt 2 ) 'rf> n ^ j > n = ^ A fa") "V&iAn • (25) 



The differential operator A (iu n ) on ipj >n leads to the dynamics of the phase angle d T 9j 
through the component, 

(o^y ■ (oj,) = ^ ■ (d^) 2 . (26) 

Hence the contribution to the dynamics of #-field due to JT coupling is given by 

5S cS = J2A (iu n ) 4*Jj,n -> ^ ■ E \ A fan) - A (0)] ■ (),.„(),. . (27) 

3,n j,n 

Now the dymanics of the orbital is determined by the propagator D(q,uj n ) defined by 

[Dfiu,*)}- 1 = U(q,u n ) + ^ • E L -^- 2 . (28) 

There are two limits of interest. In the case of weak coupling, i.e., El <C to, the dynamics 
of the orbital is determined by U(q,u n ) and the characteristic energy is of the order of to- 
Therefore we can replace uj n in Eq. (|28|) by ~ to ^ and the correction of the propagator 
is of the order of El/^. More explicitly the reduction of the kinetic energy gain AK due to 
JT coupling is estimated by replacing In II in Eq. ( p4|) by — In D in Eq. fl28|) as 

AK ~ [ du ■ I ~ E L ul/{ul + Q 2 ) , (29) 



n(g» 

where u c is the characteristic frequency of the orbital fluctuation and u c ~ to- This energy 
correction quadratically grows up with increasing u c /Q <C 1 and then saturates into the 
lattice relaxation energy El = g 2 /MQ 2 with u c /Vt ^> 1. Considering that u c ~ to ^ Q 
we conclude that AK ~ E^. As increasing E L (~ t ), we expect the saturation effect as 
AK ~ E L to/ (to + -El) as is evident from eq.(28). 

The strong couling limit, i.e., El ^> to, is more interesting. In this case, the orbital 
dynamics is determined by both II and El terms in Eq. Q2"5|) . More explicitly the propagator 
is approximated as 
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[DfauJ]- 1 = t + \ ■ E L ^f^, (30) 

and its characteristic energy scale is then given by 

Uc 9*nJlL <0. (31) 
V El 

Therefore there occurs the slow down of the orbital motion in this strong coupling limit. 
However the small polaron effect should be relevant in this case, which can not be treated in 
the present formalism, and the detailed study on this strong coupling case is left for future 
studies. 



IV. DISCUSSIONS 

We now compare the results for the non-interacting and strongly interacting limits. The 
order estimation of AKu =0 for the non-interacting limit is as follows. Rather complicated 
form of the non-interacting result, Eq. fljlf ), Eq. (|A3|) , Eq. (|A6|), and Eq. ([A 71) , roughly 
takes the ^-dependence as 



g V i , g 2 i «A 



AKu =°~ to \vm) ■ foTn) 3 • " to ' MT \%faf ~ El " (i ^ (32) 

The dependence is hence a kind of perturbative forms with the intermediate energy denom- 
inator 1/ (to + fi) 2 and the vertex g/y/MD,. The small factor f2/t comes from the fact that 
only the states with the energy window ~ Q near the Fermi energy is influenced by the el-ph 
interaction. More explicitely, 

AK AK/E L (value picked up from Fig. 4) 

~K,~ = K /E L = K /E L ' ( } 

Ei can be evaluated as ~ 0.6 eV from the literature@. With K ~ 2.16 eV in our calucula- 
tion, and the value ~ 0.1 in Fig. Q we get AK/K ~ 3% as a lower bound. 

On the other hand, in the strong correlation limit, the effect of Fermi degeneracy and 
the small factor Q/to are missing and AK ~ El- This means that the strong correlation 
enhances the JT effect. It is reported that the observed spin stiffness is well reproduced 
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even semiquantitatively by meanfleld estimations without considering the el-ph interaction 
discussed here.0IH This means that the small polaronic effect is absent in the metallic region, 
namely El is of the order or less than to (Here to should be interpreted as the "bandwidth" 
rather than the transfer integral). It is reasonable with the above estimation that El ~ 0.6 
eV. 

In summary, we have studied JT el-ph effect in three dimension with and without the 
electron correlations. In the non-interacting limit, the reduction is calculated as a function 
of the doping concentration. In this case, the doping dependence is mainly dominated by the 
density of states at the Fermi energy. It is shown that the kinetic energy is always reduced 
by the JT el-ph interaction even if the off-diagonal processes in orbital indices are taken into 
account. The reduction AK of the kinetic energy K is estimated as AK/K = 3% in this 
non- interacting case. This small value is due to the small factor £l/t occuring in the Fermi 
degenerate case. In the strong correlation limit, we have derived an effective action to study 
the el-ph interaction. The small factor fl/to is missing in this case and el-ph interaction is 
enhanced by the strong correlation. AK ~ E^t /(-Ex + to) is this case without the small 
polaron formation. From the comparison with experiments and the above results, the small 
polaron formation is unlikely in the metallic state. 

The authors would like to thank A. Millis and H. Yoshizawa for their valuable discussions. 
This work was supported by Priority Areas Grants from the Ministry of Education, Science 
and Culture of Japan. R.M. is supported by the Japan Society for the Promotion of Science 
(JSPS) for Young Scinentists during this work. 

APPENDIX A: EVALUATION OF EQ. (12) 

With afc (z, £k-q) defined by the breathing type self-energy expression^ as, 



vo ( z \ _ 9 \ - 



2 



/ (&-,) +N(u) f + N {-uj) 



the JT type self-energy Eq. (11) can be written as 



E > ( A1 ) 
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S F' (*) =E [Onl-g + 4S-g) 0> S l-J + 0^1-* + ^J-J (*> S U)] > (A2) 

with coefficients defined in Eq. (§). N (u) represents the bose distribution function. 
The notation ^ ±;n) is defined as £ k ±;n) = + Q. Substituting Eq. @ into Eq. (pD 
leads to 

j (ti70 

(727^7370 _ p 7i7i(+) 4 727^ 4 737^ X T, (-(+) 

^jfc-g ' v 5 " k 1 ^k-q 

, P 7i7i(+) 4 727a 4 7373 vT/ /->(-) ~(-) ~(- 



g 2 /2NMu 



(+) 



, ( ei717i(+) 4 727a 4 737 3 , P 7i7i(+) 4727 2 4 737 3 \ xT/ /->(+) ~H ~(+) 
K^k-q ^-jfe ^fc-g A -\h /i +;k J ' W \^k > "fe > "fc-g 

, P 7i7i(-) 4 727 2 4 737 3 vT/ /->(+) ~(+) -•(-) 



+;fc +;fe ' 1—1 ^ ' ^k—q 

, P 7i7i(-) 4 727 2 4 737 3 , T/ /->(-) ~(-) ~H 

, (ffnUli-) 47272 4 737 3 , R 717U-) 4727 2 4 737 3 \ . ^ ~(-)\ /ao\ 

where we defined 

Ek-q^ ^ = ^±;fe-g + ^±;fe-g 5 (A4) 

[the first (second) term corresponds to the scattering by u- (v-) phonon, respectively, as in 
Fig. 2]. Function \1/ is defined as the integral 



corresponding to a diagram shown in Fig. |5|. 



T > hk-q J .... 

M 1 ) 




FIG. 5. Diagram corresponding to * , ^ 2 "* , C^-q) 
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Depending on the degree of the pole, it is evaluated as, 



(3) 



/ 




/ 
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(^) 








(ff-ff) 


(<?> 





/(<er') 



(«Er" -<£")(« 

/(«. 



(3;-u) _ ,(2) 



r(i) 



+ 



Ml) _ t(3;+a/)\ / Al) _ A2)\ 1 /.(2) _ (A2) _ A3;+u>)\ 

\<ik Sfc- 9 J \$k Sfc ) \^k Sfc y l^Sfc Sfc-g J 



f (£r } ) 



/t(3;+^) _ Ai)\ (A3;+u) _ A2) 

[yk-q Sfc J [yk-q ^k 



(A6) 



for + if and 



>k—q 



f (4) (ik - - / (4) + / 



t e+u 



+/(& 



^ /' (£fc) Cfe-g, 


-/(&)+/ 








- ik- qj 


2 



-/ (6k) 



fUk-q fUk-q IP 



(Zk £k-q) 
+ f(dk~q)f(ti-q)-P 



+/(&) 



/ (ik-q) f (C-J ■ P 



(4 - e 



—a; 



(4 - e 



A;— 5 



(1,2) . , T .(3,4) 



(A7) 



where /(£fc) — 1 — f With notations here, the expression for the simple breathing type 
case with the single band systemll is expressed as 

,2 



ANMu 



(A8) 



kq 
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APPENDIX B: PROPERTY OF THE FUNCTION * IN EQ. (A7) 



For the simplest case with single-band electrons and breathing type phonons, Eq. (KE) 



the filling dependence is determined by [~^k 4 ^ (6fcj 6fc> Each term of the function 

kq 



^ £jfcj £fc-g) m Eq. (|A7|) contributes when the states |^), |£fc- 9 ), and |^_ 9 ±a;) are in 



such configurations as shown in Fig. ^](a). 



If.- 



1st & 2nd terms; k e occupied 



1^} 



1st term 



If.) 



2nd term 




3rd & 4th terms; k e unoccupied 



14) 



If..) 



I?.) 



3rd term 



4th term 




(b) 



FIG. 6. Configurations of |£fc), |£fc- g ), and under which each term of Eq. ( A7) 

contributes. 



For given and fixed fc, the contribution due to each term behaves as shown in Fig. |7] as a 

( 1 2^ \ 

function of Notice that the first and the second terms (W k ' ) are exclusive to the third 
and fourth terms (^i 3 ^) each other. Introducing a notation k = k — 2k f [see Fig. ||] (b)], 
Eq. (|A8|) is evaluated as, 

fcg fcg k£occupied V 9 <? 

= E (e fc -^ 2) +e,-^ 4) ), (Bl) 

kdoccupied 

where > and Sf A) < are the quadratures of the shaded with signs depicted 
in Fig. [7| (a) and (b), respectively. 
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k- q 



-CO 
1st term 



- ■ - 



T 



ferm 



0) 



- 6., 



3rd term 



\ 



A 



4-, 



ft:- q 

(P) 



FIG. 7. Behavior of each contribution of terms in Eq. (A7) as a function of £,k-q for fixed k. 



(3,4) 



(1 2) 

SI + 5Sk with a small deviation 5Sk which reflects the 



It can be written as S 
difference between the band curvature at \k) and k^J shown in Fig. |B| (b). 5Sk vanishes when 
the system is half-filled where the fermi level locates at the middle of the band. Equation 
rTT) is then evaluated as, 



AK sin 9 le ^J2[(H + &) S? A + (/i - &) ("4 1,2) + SS„ 
kGoccupied 



kGoccupied 



(B2) 



with a small 5K compared with the first term. Because < for occupied states the result 
has therefore definite sign. From Fig. [7|, it is understood that the contributions mainly come 
from the vicinity of the fermi level. It thus means that the result scales to the density of 
states at the fermi level. 
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